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Abstract 



In this paper we compute, within in the context of a relativistic quark 
model, the Isgur-Wise functions for exclusive semileptonic B Xc decays, 
where Xc is any charmed mesons with total spin J = 0, 1,2 or one of their 
first excited states. The relevant matrix elements are computed by a direct 
numerical integration, in coordinate space, of the convolution of the wave 
function of the B meson at rest and the wave function of the Xc meson, 
boosted according with its recoil factor. Our results are compared with other 
predictions found in the existing literature. 
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I. INTRODUCTION 



There has been much progress in recent years in model-independent calculations of heavy 
meson decays. By using Heavy Quark Effective Theory (HQET) [0] and Lattice Gauge The- 
ory, we can now make some definite predictions for certain processes, in a limited kinematic 
range. Knowledge of these decays are extremely important for particle physics, not just in 
their own right, but in measuring fundamental parameters such as Vcb and sin(2/5). 

Unfortunately, at present time, it is not always possible to use these techniques for all 
kinematic situations. One example are the Isgur-Wise (IW) functions, which relate all the 
different form factors for heavy-to-heavy decays to a single function, at leading order in 
A/m. These functions can, in principle, be determined by lattice computations but, because 
of limited computing resources, they have only been computed for decays into the ground 
state and the error is still sizable. For the moment HQET can be used to produce model 
independent results at zero recoil but, away from this point, the IW functions are unknown. 
Therefore, for the time being, we must rely on models. 

We present here a study of B decaying into excited Ds and a determination of the 
corresponding IW functions, following the work of Ref. P|. Our study is based on the quark 
model proposed in Ref. [Q], where a Dirac equation was used to describe the light quark in 
the potential of the heavy quark and determine masses and wave functions of excited mesons. 
We use these wavefunctions to calculate the leading order IW functions for B Xc decays, 
where Xc is a spin 0-2 charmed meson or its first radially excited state. The IW functions 
are computed explicitly by a three dimensional numerical integration of the relevant matrix 
elements expressed in terms of the wave functions derived in Ref. 

The paper is organized as follows. In Section II we set up the formalism and in Section 
HI we discuss the actual calculation of the IW functions. In Section IV we discuss the results 
and compare with the literature. Finally in Section V we conclude. 
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• P TT; 

J = Si 


Particles 




m (GeV) 


m! (GeV) 


1- 
2 


D 


0" 


1.865 


(2.589) 


1 + 
2 


D* 


1" 


2.007 


(2.692) 


n* 


0+ 


(2.377) 


(2.949) 




Di 




(2.490) 


(3.045) 


3 + 
2 


Di 


1+ 


2.422 


(2.995) 




Dl 


2+ 


2.459 


(3.035) 


3- 
2 


Dl* 


1" 


(2.795) 


(3.420) 




D*2* 


2" 


(2.833) 


(3.459) 



TABLE I. Charmed meson spin multiplets (g = u, d). The masses are experimental when 
possible, otherwise they are calculated from the model. The primed masses are for the first radial 
excited states. The masses in parenthesis are predictions of the in the model. 



II. THE QUARK MODEL 

Due to heavy quark symmetry, we know a great deal about the spectroscopy of mesons 
containing a heavy quark Q. At leading order in K/itlq, the spin and parity of the light 
degrees of freedom, ' , are good quantum numbers. Thus the particles appear in multiplets 
labeled by both the total spin J and In this limit, the mesons come in degenerate 
doublets with total spin Jj. = s; ± sg, where sq = 1/2 is the heavy quark spin, a conserved 
quantum number at leading order in HQET. In Table | we list the low spin charmed mesons 
considered in this paper, their masses and corresponding quantum numbers in HQET. In 
those cases where the experimental mass is not known we report, in brackets, the mass 
predicted in Ref. For the mass of the decaying B meson we adopt the PDG value 
TUB = 5.279 GeV. 

The model of Ref. is a relativistic quark model characterized by a spin-dependent 
plus a spin-independent potential. The general form for these potentials has been derived 
by general arguments such as confinement and asymptotic freedom. The parameters that 
appear in the potential have been fixed by fitting the experimental spectrum of Z^, Ds, B and 
Bg mesons with the predicted spectrum, including 1 / mf, corrections |^ and mixing effects. 

The most general heavy-light meson (in the D^, B, Bg family), H, is a bound state of 
a light quark (g) and a heavy quark (Q). The heavy quark is treated as a static source of 
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chromoelectric field and the only quantum number associated with it is its spin. The light 
quark is treated relativistically and its state is described by the wavefunction ipn,e,j,m{.r, 9, ip). 
We introduce the following quantum numbers: 

• n, the number associated with the radial excitations 

• the orbital angular momentum 

• j, a short hand notation for S£, the total angular momentum of the light quark 

• m, the component of j along the z axis 

• J, the total angular momentum of the system 

• M, the component of J along the z axis 

• 5, the spin of the heavy quark along the z axis 

The parameters of the model are the masses of the light quarks (m^ for q = u,d or s), the 
masses of the heavy quarks {rriQ for Q = c or b) and the chromoelectric potential of the 
heavy quark {V{r)). 

The total wavefunction of the system can be decomposed as follows 



(1) 



where C^'^i „ are the usual Clebsh-Gordan coefficients and fs is a two component spinor 
representing the heavy quark. In Eq. ([l|), the four spin components of the light quark 
wavefunction, is parametrized as follows: 



o . 1 „ 



2 

1 

2 



(2) 



Here Y^{9,ip) are spherical harmonics that encode the angular dependence while /° 
fn^ijir) are real functions that encode that radial dependence, kf^j^m ^'^^d k^^^^ are fixed, up 
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as 


A 


b 


m„ 


rus 






TTlb 


Mb 


0.339 


2.823 


0.257 


0.071 


0.216 


1.511 


1.292 


4.655 


4.685 



TABLE II. Parameters for the model. The masses and A are all measured in GeV. The param- 
eter b is measured in GeV'^. 

to an overall phase, by imposing a normalization condition. Our choice of the phase is such 
that 



The Hamiltonian of the most general heavy-light system, at leading order in l/rrih reads 

=7°(-z^ + m,) + F(r), (4) 

and the rotational-invariant potential is the sum of a constant factor (Mg), a scalar part 
{Vs) and (the zeroth component of) a vector part (l^) 

V^(r) = MQ + 7°K(r) + V;(r), (5) 

where 

4 ry 

K(r) = --^erf(Ar), (6) 
Vs{r) = br + c. (7) 

The role of the erf() in the potential is that of regularizing the ultraviolet divergence in the 
l/rrib corrections to the spectrum. The parameters of the model, determined in Ref. [Q, are 
shown in Table |I|. 

III. CALCULATION 

The hadronic part of the semileptonic exclusive decay B D + i + u (for the 
most general excited D in the final state) is encoded in a matrix element of the form 



{D{n',i',j',J',M')\r\B{n,£,j,J,M)) where f = Uc{0)rub{0). For the decays of interest 
r = or 7'^7^ but for the purpose of this paper we consider the most general F structure. 
In fact, thanks to the HQET, all matrix elements that differ only for the spin structure can 
be related to the same IW function. The heavy-light states are normalized according with 
the usual non-relativistic convention. 

A general formalism for the computation of matrix elements between states represented 
by wave functions was derived in Ref. 0. In that paper the problem of defining equal time 
wave functions is discussed and the matrix elements are written as a integral in momentum 
space of the Fourier transformed wave functions. The general problem is greatly simplified 
in the specific case of heavy-light systems since it is always possible to shift a meson in time 
by changing the phase of the heavy quark. Hence, in this paper, we find more intuitive to 
express our matrix elements as integrals in coordinate space which, in the most general case 
of interest, look like: 

{D{n',e,f,J',M')\T\B{n,i,j,J,M)) = J ^illl,^^,^j,^,,,{^)t<iI^,e,,,j,Mi^) d'^, (8) 
where is the wavefunction \I' boosted according with the recoil momentum p: 

Ke.,j,M(^) = E Cfl'^ s<e,,mi^) ® ^s. (9) 

and 

^P(x) = 5(Ap)^(Apix), (10) 

es = S{A^)^s. (11) 

where Ap is a boost in direction p. 

Note for any state only that component of the spin parallel to the direction of motion 
p is a good conserved quantum number, the helicity. Therefore in Eq. we chose m', M' 
and S to be the components of the light angular momentum, the total angular momentum 
and the heavy quark spin respectively, parallel to the direction of the boost Ap. We checked 
that, with this definition, the result for the matrix element is independent on the direction 
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of p. In our analysis we ignore mixing in the wavef unctions and other 0(m^^) corrections.^ 
The 3D integrals are evaluated numerically using the Vegas Monte Carlo algorithm. 

In order to compare matrix elements we calculate with the model to the corresponding 
IW functions, we need to calculate the matrix elements in HQET. This can be done by using 
the trace formalism As an example, consider the s^' = | doublet. The fields Py and 

P*^ that destroy the members of this doublet with four-velocity v are grouped together in 
the 4x4 matrix 



H- 



1+' 



(12) 



The matrix satisfies the relations 'f)H^ = = —H^ f. To leading order in Aqco/'^b.c 
and OLsi the matrix element between B^*^ and Z^^*-* mesons are 



(13) 



Here ^(w = f ■ f') is the (dimensionless) IW function for B decaying to the D*^*^ multiplet. 
Matrix elements with any F can now easily be calculated and are related by heavy quark 
symmetry to £^{w). 

For other multiplets, different 4x4 matrices are used. The general form for arbitrary 
spin was developed in Ref . |^ . Below are the other matrices which are necessary to calculate 
relations used in this paper: 



V 



1 + 



p: '^tmt' - p. 



(14) 
(15) 
(16) 



^Ref. ||5| suggests that l/m;, corrections to the matrix element corresponding to ^ D^*^ are 
smaller than expected by dimensional analysis. We do not know if this is also the case for the 
other matrix elements of relevance in this paper. 
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FIG. 1. IW function for the D — D* multiplet, as predicted by different models. The vertical 
lines mark the end of the kinematical allowed region of the two mesons in the doublet. 

where is for the sf' = |^ doublet, and F^'^ are for the sf' = |^ doublets. 

Using the trace formalism, we can relate the matrix elements calculated in the model, 
Eq. (P), to the IW functions. Due to heavy quark symmetry, there are many matrix elements 
that could be used to obtain the same IW function. By using different choices, we can check 
to make sure the model is giving consistent results. Also, many matrix elements are equal 
to zero at leading order in AQCD/''^fe,c and a^, which is another way to check the model 
results. In the Appendix, we show the relevant matrix elements for the different doublets 
and different spin structure. As an example, again consider the sf' = \ doublet. Picking 
r = 1 we have 

e(^) = TT^^^I^^^I^) = TT^ / ^1, 0. i 0. W 1 ^1, 0, 1 0, o(x), (17) 
and similar relations can be found for other doublets. 

IV. RESULTS 
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A. The D and D* multiplet 



The most investigated heavy-to-heavy decays are B ^ D and B ^ D* which corresponds 
to the s'l'- = I doublet and are parametrized in terms of the same IW function ^{w): 

aw AStt-^ 

y w + 1 [mB + niD*) J 
We know from heavy quark symmetry that, ignoring perturbative corrections, the IW func- 
tion for these decays, ^, is normahzed to one at zero recoiL \Vcb\ can be determined from 
exclusive B D^*^^ decay channels by a direct application of Eqs. ( [T^JT^ ). On the exper- 
imental side it is difficult to extrapolate C,{w) at zero recoil, since the corresponding matrix 
element vanishes as (w^ — 1)^/^ for D* and (w^ — 1)'^/^ for D. However, analyticity imposes 
stringent constraints P,p!0[|. 

The normalization of C, at zero recoil is readily visible in Fig. |l], where we plot ^ as a 
function of w = v ■ v'. If we parameterized the shape of ^ as 

e(l) [l - p\w - 1) + c{w -ly + ...], (20) 

we obtain for the slope parameter = 0.501 and for the curvature c = 0.145.0 In Fig. |l| 
the same IW function as predicted by the following alternative models: 

• ISGW [0]. In this model, nonrelativistic meson wave functions are obtained using a 
variational approach to the Schroedinger problem and approximated with harmonic- 
oscillator wave functions. 



SBW ||TJ. In this model, the form factors are calculated assuming a pole structure. 



Note that our result is not in agreement with Uraltsev's sum rule p > 3/4 
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Isgur-Wise function 

B(n=1, L=0, j=1/2) -> D(n=1, U1, j=3/2) 



0.14 



0.12 



0.10 - 



0.08 




06 - ■ ■ Model 



0.04 



FIG. 2. IW function for the Di — multiplet. The vertical hnes mark the end of the kine- 
matical allowed region of the two mesons in the doublet. 



Covariant Oscillator Quark Model or COQM This model is based on a covariant 
representation for nonrelativistic meson wavefunctions. 



If Fig. ID we also compare our model with the lattice prediction of Ref. |jT5|. The lattice result 
is affected by unknown systematic quenching errors and discretization errors, particularly 
for large momentum transferred. These errors are difficult to estimate at present and are 
not reported in our plot.0 



We also observe that a direct lattice determinations of the slope this IW exists |jT6|. This 
computation is done with a propagating heavy quark slightly heavier than a charm meson 
and a light quark with a mass about the strange mass. It predicts a value of = 1.7 ±0.02, 
which is about a factor two larger than quark model predictions. 



10 



B. The Di and Dl multiplet 

The narrow resonances Di — D2 with hght quantum numbers sj^ = |^ are important 
for a number of reasons. For example, it is interesting to understand the composition of 
the inclusive B semileptonic decay rate in terms of exclusive final states. The particles in 
the 1^ doublet are important exclusive channels for this comparison. It is also important to 
know the decay spectrum for these particles as one the dominant backgrounds to i? ^ D'-*^ 
decays. Finally, there has been renewed effort in constraining the slope parameter for 
B D^*^ using sum rules and data on B decays to excited D states [r7,rT|. 



In Fig. I we plot the IW function, r3/2, for the D1—D2 multiplet. Unlike the previous case, 
there is no reason why r3/2 should be normalized to one at zero recoil. If we parameterized 
the shape of r3/2 as t^/2{w) = r3/2(l)[l - P3/2(^ - 1) + C3/2{w - 1)^ + ■■•], we obtain 
for the normalization r3/2(l) = 0.122, for the slope parameter pg^g = 1.171, and for the 
curvature C3/2 = 0.601. In Table |IT| we compare our result with other predictions found in 

P3/2 

r3/2(l) is lower. 



the literature. Note that our our value of Pn,^ is consistent with the other models, while our 



C. The and Dl multiplet 

Only recently data have been produced on the rates to these excited mesons The 



sj'' = I multiplet is very broad, as it can decay strongly to D^*'^n in an S wave |^ (the 
current experimental width is 290^79^ ± 26 ± 36 keV for the J = 1 meson), while Di and D2 
can only decay through a D wave, thus being narrower resonances (the experimental width 
is lS.9tti keV for the J = 1 meson) p^. 

In Fig. ^ we plot the IW function, T1/2, for the sj^ = Dq — Dl multiplet. Again, there 
is no reason why T1/2 should be normalized to one at zero recoil. If we parameterize the shape 



^Modern computer technology allows for a noticeable improvement over these lattice results. 
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Ref 



Ours 



^3/2(1) 
0.12 
0.49 

0.41 
0.56 
0.66 

0.54 
0.52 
0.31 



^3/2 
1.17 

1.53 

1.5 

2.3 

1.9 

1.5 
1.45 
2.8 



Tl/2(1) 

0.094 
0.28 
0.41 
0.09 
0.41 
0.35 ± 0.08 
0.22 
0.06 
0.31 



Pl/2 

0.821 
1.04 

1.0 

1.1 

1.4 
2.5 ± 1.0 
0.83 
0.73 

2.8 



TABLE III. Comparison of IW functions r3/2 and at zero recoil and their respective slopes 
P3/2 and p2^2 



and p?/o from different models. 



of ri/2 as ri/2(w) = ti/2(1)[1 — Pi/2('^~l) + Ci/2(^ — 1)^ + - • we obtain for the normaUzation 
^1/2(1) = 0.094, for the slope parameter p^j^ = 0.821, and for the curvature ci/2 = 0.244. In 
Table |IT|, we compare our result with other predictions found in the literature. 



D. The Dl* and multiplet 

To be complete, here we briefly mention the only other doublet containing a spin one 
meson, the = | , corresponding to L = 2 orbital excitations in the quark model. The 
states are expected to be even broader than the sf' = ^'^ multiplet. In Fig. || we show 
the IW function, k for this doublet. If we parameterize the shape of k{w) as k,{w) = 
k(1)[1 — p^fli'^ 1) + C3/2('W^ — 1)^ + . . .], we obtain for the normalization = 0.367, for 
the slope pyl = 0.884, and for the curvature = 0.377. 
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0.100 
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0.060 



0.040 

1 1.1 1.2 1.3 1.4 

w 

FIG. 3. IW function for the — D\ multiplet. The vertical hnes mark the end of the kine- 
matical allowed region of the two mesons in the doublet. 

E. The Radial Excitations 

While none have been seen, there are radial excitations of all the above mentioned 
doublets. It is unlikely that they will be seen in the near future, but their effects are 
important in reconciling the inclusive b —* c and exclusive B ^ Xc semileptonic decay 
rates. The radial excitations are also important as they enter into sum rule calculations. 
We therefore discuss the first radial excitations of the above doublets here. We denote the 
IW functions for the radial excitations with primed versions of the same Greek symbols as 
their non-radially excited counterparts. 

At zero recoil, the IW function for B D^*^ was normalized, ^(1) = 1. For the radial 
excitation, there is no overlap between the B and D'^*'^ at zero recoil, thus ^'(1) = 0. This 
can be seen in Fig. |^. If we expand ^' around w — 1, we get 

^'(w) = -0.325(w - 1) + 0.213(w - if. (21) 

Our result is very different from the conclusion of Ref. They get C,'{w) = 2.2{w — 1) + 
2.6(w- 1)2. 

For the radial excitations of the other doublets sj' = |^ we do not have any nor- 
malization requirements at zero recoil (since the matrix elements vanish at leading order 
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0.40 




0.25 I ' ' ' ' ' 1 

1 1.1 1.2 1.3 



FIG. 4. IW function for the DJ* - D2* multiplet. The vertical hnes mark the end of the 
kinematical allowed region of the two mesons in the doublet. 

because of the heavy quark spin symmetry). In Fig. |^ we show the IW functions for 
t[^2 respectively. 



V. CONCLUSION 

Information on the decays of B mesons into charmed mesons is important for a number 
of reasons. Using Heavy Quark Effective Theory, these decays can be written in terms 
of Isgur-Wise functions which parameterized the form factors. These functions cannot be 
calculated except, eventually, on the lattice. 

In this paper we derived the Isgur-Wise functions of B ^ Xc within a quark model, 
where Xc can be any spin 0-2 charmed meson or one of their first radially excited states. 
Our results are compared with independent predictions found in the literature. Our model 
differs from the others because we derived the wavefunctions of the B meson and excited 
charmed mesons in a relativistic fashion by fitting the experimental spectrum with model 
predictions. 

Our results for the Isgur-Wise function, ^, for B to D^*^ decays is consistent with other 
model predictions but its slope is milder. Our result is also consistent with the prelimi- 
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FIG. 5. IW function for the first radial excitations with spin ,1\=2. The vertical lines mark the 
end of the kinematical allowed region of the two mesons in each doublet. 



nary lattice results of Ref. ||T5[. For decays into the P- waves the situation is even more 
uncertain since different models predict a wide varieties of results. We believe that further 
model independent investigations are required in order to put more constraints important 
phenomenological quantities. 
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APPENDIX A: MATRIX ELEMENTS IN HQET 



In this Appendix we collect the relations between the matrix elements of the form 
{D\cTd\B) to expressions in terms of the Isgur-Wise functions in HQET at leading or- 
der in AQCD/''^f),c and a^, in analogy with Eq. (0). The matrices containing the states for 
the different doublets are shown in Eqs. ([1^) and ([l^p!6D . The choices of Dirac structure 
are F = (1, 75, 7^^, 7^75). In this appendix is the polarization vector for spin 1 particles 
and e^^ is the polarization tensor for spin 2 particles, while e"^"'^ is the usual Minkowskian 
antisymmetric tensor, v is the velocity of the decaying B meson and v' is the velocity of the 
charmed decay product. All matrix elements not explicitly shown are zero: 



{ D{l/2-) I cYb 

{ D*{l/2-) I ci'h 
{D*{l/2-)\ci^b 
{D*{l/2-) I cYl'b 

{ /^o(V2+) I n'b 
( I^o(V2+) I c^l'b 

(A*(V2+)|c6 
(Dt(l/2+) |C7^6 
(Dt(l/2+) IC7V& 

( Di(3/2+) I cb 
(A(3/2+) |C7'^6 



B 
B 

B 
B 
B 

B 
B 

B 
B 
B 

B 
B 



{l + w)^{w) 

(v+v'niw) 

[e^\l + w)-e''v,v'^]^{w) 
2{w - l)ri/2(w) 

-2e''vuTi/2{w) 

2[e''(u^-l)-eW1ri/2H 



-y2eX(l + w)r3/2(ti;) 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 

(A6) 
(A7) 

(A8) 
(A9) 
(AlO) 

(All) 



^[£^^(1 - w^) - 'it^'v.v^ + (w; - 2Yv,v'^ T^i2{w) (A12) 
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( I)i(3/2+) I n'l'b I 5 ) = --^(1 + w)e^^P^e^vpv'^n/2{w) (A13) 

( /?2*(V2+) \ci''b\B) = ^/3e"V^;;3T-3/2(u;) (A14) 

( L';(3/2+) I C7'^M ^ ) = -iVSe'^'^P'^eaaVvpv'^nMw) (A15) 

( /?2*(3/2+) I I ^ ) = i;,[e"^t;^^;"' - v^(l + w)e"^]r3/2(w^) (A16) 

( Di(3/2-) \c^%\B) = y|(l - «;)e^^;,fi;(w;) (A17) 

( D,(y2-) I C7^M ^ ) = - w)e''''^^e^Vfsv'^K(w) (A18) 

( L'llS/s") I cYl'b \B)^ /|^1(2 + ^^)^-^^''' - S^-^'' + (1 - ^')^^^]'«(«^) (A19) 

( L'2*(3/2-) I cM ^ ) = ^;«^;;3e"^/^(«;) (A20) 

( ^2(^2") I C7^6 I ) = [(1 - w)e"^v^ + e'''^v'^Vo,vp]K{w) (A21) 

( D;{3/2-) I C7'^7'^' I ^ ) = ^e'^"^^e«,^XAc(^) (A22) 
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